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This report looks at the interaction of radiated electromagnetic fields with 
earth ground in military or law-enforcement applications of high-power 
microwave (HPM) systems. For such systems to be effective, the microwave 
power density on target must be maximized. The destructive and construc- 
tive scattering of the fields as they propagate to the target will determine 
the power density at the target for a given source. The question of field 
polarization arises in designing an antenna for an HPM system. Should 
the transmitting antenna produce vertically, horizontally, or circularly po- 
larized fields? Which polarization maximizes the power density on target? 
This report provides a partial answer to these questions. The problems of 
calculating the reflection of uniform plane wave fields from a homogeneous 
boundary and calculating the fields from a finite source local to a perfectly 
conducting boundary are relatively straightforward. However, when the 
source is local to a general homogeneous plane boundary, the solution can- 
not be expressed in closed form. An approximation usually of the form 
of an asymptotic expansion results. Calculations of the fields are provided 
for various source and target locations for the frequencies of interest. The 
conclusion is drawn that the resultant vertical field from an appropriately 
oriented source antenna located near and above the ground can be signif- 
icantly larger than a horizontally polarized field radiated from the same 
location at a 1.3 GHz frequency at observer locations near and above the 
ground. 
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1. Introduction 


Effective military or law-enforcement applications of high-power mi- 
crowave (HPM) systems in which the HPM system and the target system 
are on or near the ground or water require that the microwave power den- 
sity on target be maximized. The power density at the target for a given 
source will depend on the destructive and constructive scattering of the 
fields as they propagate to the target. Antenna design for an HPM sys- 
tem includes addressing the following questions about field polarization: 
Should the fields the transmitting antenna produces be vertically, hori- 
zontally, or circularly polarized? Which polarization maximizes the power 
density on target? (The question of which polarization best couples to the 
target is beyond the scope of this report.) While this report does not com- 
pletely answer these questions, it addresses the interaction of the radiated 
electromagnetic fields with earth ground. It is assumed that the transmit- 
ting antenna and the target (or receiver) are located above, but near the 
surface of a flat idealized earth (constant permittivity, e, and conductivity, 
cr) ground. First an ideal vertical dipole (oriented along the z-axis perpen- 
dicular to the ground plane) is addressed. The horizontal dipole (parallel 
to the ground plane) follows. 
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2. Problem Formulation 


The problems of calculating the reflection of uniform plane wave fields 
from a homogeneous boundary and calculating the fields from a finite 
source local to a perfectly conducting boundary are relatively straightfor- 
ward. However, when the source is local to a general homogeneous plane 
boundary, it is found that the solution cannot be expressed in closed form. 
An approximation usually of the form of an asymptotic expansion results. 
The problem of an ideal dipole over a homogeneous half-space has been the 
topic of a number of studies starting at the turn of the last century with the 
solution provided by Sommerfeld (1949) and leading to the more contem- 
porary work of Banos (1966) and King et al (1994, 1992). References such 
as Maclean and Wu (1993) address in detail the many approaches to solv- 
ing the problem. Considerable controversy has surrounded these studies. 
We will not attempt to derive the solution or otherwise discuss the solu- 
tion of the problem in this report. We will rely on the work of King et al 
for a complete and concise formulation of the problem. Figure 1 depicts the 
problem geometry. The King expressions have been encoded and solved in 
MATLAB® (1984-1999). Comparisons of the field structures are provided 
for various source and target locations for frequencies of interest. 

The expressions that follow are constrained by the magnitude of the wave 
numbers (k = ta(pe) 1//2 , ui = 2ir x frequency) in each region: 

M>3|fc 2 |, (1) 

where, after King, the subscript 2 denotes the upper half-space (air) and the 
subscript 1 denotes the lower half-space (earth). Also, // is the permeability 
in henries per meter, e is the permittivity in farads per meter, and a is the 
conductivity in siemens per meter. The subscripts 0 and r represent free 
space and relative to free space, respectively. 
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2.1 Vertical Dipole Over Earth 

The electromagnetic fields from a dipole with dipole moment II oriented 
perpendicular to and at a height z = d above the ground plane have three 
components in cylindrical coordinates. The magnetic field is symmetric 
about the z-axis, perpendicular to the direction of propagation. For that 
reason the fields will be termed transverse magnetic (TM) fields. (We will 
find that the horizontal dipole has TM and transverse electric (TE) com- 
ponents.) The formulation as provided in King et al (1994) is, referring to 
figure 1, 
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and C 2 (P), S 2 ( P ) are the Fresnel integrals as defined in Abramowitz and 
Stegun (1970). Since 


C 2 (P) + iS 2 (P)= 1 -(l + i) 



( 8 ) 


we have 



(9) 


where w(z), called the plasma dispersion function with complex argument, 
is a form of the error function defined in Abramowitz and Stegun (1970). A 
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convenient MATLAB® m-file is available for solving the function with com- 
plex argument (Chase) and is provided in the appendix. Note that \fi can 
be written as ( ? 1 ^ and that F(P) always appears with e~ lP , thus cancel- 
ing the exponential term in equation (9). This is reflected in the MATLAB® 
m-files that calculate the fields (provided in the appendix). 

The equations have been written such that the direct and reflected compo- 
nents appear first. The last term is referred to as the surface or lateral wave. 
Often it is called the Norton surface wave from the engineering models he 
developed in the mid-1930s. The equations reduce to the fields above a per- 
fectly conducting ground when k\ — > oo; the surface or lateral wave term 
then goes to zero. 

2.2 Horizontal Dipole Over Earth 

The electromagnetic fields produced by an ideal dipole oriented at a height 
z = d above and parallel to the ground plane consist in general of both TM 
and TE components. The formulation as provided in King et al (1992) for 
the TM wave components in cylindrical coordinates, referring to figure 1 
and the above definitions, is 




U A 
— cos 0 

47T 


Jk 2 n ( z - 


( z~d \ f iki _ 1 1 _ „ik 2 r 2 ( z+tA f ik 2 _ 1 
V r i ) V r i if/ ' \ r 2 ) \r 2 r| 


2 k 2 „ik 2 ro. 
' fcl e 


a ik 2 r i 


ik 2 
r i 
ik 2 
r 2 


1 

pr 

'2 


k 2 rS 


F p {p,4>,z) = 


ojHoH 

4irk2 


COS 0 < 


— e 


ik 2 r 2 


[-(l) (?) fe)^- ip F(P)_ 

( z~d \ ^ ( OSZ 

\ ri ) \ri 

m 2 ( 

(^)( 


2 , 2 i , 

if k 2 r\ 


2 2 i 

if ^ *fif 

2 k 2 
ki 


3 i 
k 2 r 


ik 2 

r 2 


3 

pr 

'2 


l) 

Si 
k 2 r ^ 


| 


ik 2 

r 2 


ik 2 

r 2 

1 

pr ■ 

'2 


1 

pr 

'2 


k 2 r% 


_ f M') ( T2.\ ( ^ 2 e —iP 


E z (p,(j>,z) = 


u>po n 
4irk2 


cos 0 < 


-e ibn (if) (t r) (ff 
+e iter > (if) (f? ? 


( ^2 ^ ( V2. ^ ( 

\ki J \ P J \ k2V2 J 

?) 


F(P) 


F2 


3 

3 

pr 

'2 


Si 

k 2 r 

Si 

k 2 r% 


2fc 2 „ik 2 r 2 
kl e 


(if) (t - A) - § (i^N- iP np) 


) 

7T 

k 2 r 2 J 


(10) 


> ,and (11) 


> • 


( 12 ) 


The TM fields are zero broadside to the dipole orientation, d = f • The TE 
components are 
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Hp(p, (j), Z ) 


H z {p, 4>i z ) 


E<p(p, Z ) 



These are the predominant fields broadside to the dipole orientation. Here 
Eff, is often termed the horizontal electric field. Again, the equations reduce 
to the fields above a perfectly conducting ground when k\ — oo; the surface 
or lateral wave term goes to zero. 


2.3 Comparisons 


The peak power radiated by a unit dipole is given in Collin and Zucker 
(1969): 


p fcjCo 
"127 r ’ 


(16) 


where Co is the free-space impedance ■ The fields and power density 

comparisons that follow are normalized to one watt radiated peak power. 
To obtain the unnormalized quantities, simply multiply the power results 
by P and the field results by Pz. The outward component of the complex 
Poynting vector over a closed surface is 



■dS = - P 


complex 1 
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where the negative sign indicates power flow away from the surface. Our 
interest is in the real part of the power density at the observer (or target) 
location. In our cylindrical coordinate system for the TM components, this 
becomes 


V2 (E x H* ) = 1/2 ^ (p^H^ * + z 0 E p H^ *) , (18) 
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where p and z 0 are the unit vectors in cylindrical coordinates. Near the 
ground, the power flow is predominantly radial with a small z component. 
And, for the TE components, we have 

v 2 Re (E x H*) = 1/2 Re (p^H, * - z 0 E^H p *) . (19) 

For the vertical dipole, the power density at the observer (on target) will be 

Pv = V2 (Re {E t f,H z *)+ Re {E^H p * ) ) . (20) 

The horizontal dipole in general will produce a target power density of 

P h = I /2 ({Re (E 4> H Z *) - Re (E Z H 4> *)} + {Re (E p H 4> *) - Re (E 4> H p *)}) . (21) 

For broadside calculations this becomes 

Ph = V 2 (Re (E 4> H z *) + Re (E 4) H p * ) ) . (22) 

The calculations that follow are limited to a frequency of 1.3 GHz (wave- 
length, Ao, is 0.23 m) and dipole heights of 1 to 3 m. Observer (target) 
heights range from 0 to 5 m. The frequency of 1.3 GHz is chosen, since most 
of the HPM source and antenna design work at ART is centered around that 
frequency. The height ranges are chosen to be consistent with ground vehi- 
cle source and target applications. Five classes of ground parameters that 
are representative of distinctly different terrain will be addressed. These 
five classes are those discussed in King et al (1994) and are given in table 1. 


2.3.1 Comparison of Field Components Near the Earth 

When thinking about the interaction of electromagnetic waves with the 
earth's surface, we often view the earth as a perfect conductor. The hori- 
zontally polarized field is reflected with the opposite sign of the incident 
field, leading to a difference or destructively interfered-with field. The ver- 
tically polarized field is reflected with the same sign as the incident field, 
leading to a sum or constructively interfered-with field. Figure 2 compares 
the resultant primary field components produced from a horizontal and a 
vertical ideal dipole located at the same point in space over a perfectly con- 
ducting ground plane. From these calculations, one might conclude that 


Table 1. Earth 
parameters. 


Case 


a* S/m 

Er* 

1 

Sea water 

4 

80 

2 

Wet earth 

.4 

12 

3 

Dry earth 

.04 

8 

4 

Lake water 

.004 

80 

5 

Dry sand 

.000 

2 


*The variable name representing a 
in the MATLAB m-files is SIGMA 
and for e r it is EPSREL. These vari- 
able names appear on many of the 
figures. 


6 



Figure 2. Comparison of 3 
main electric field 
components from a 
vertical ( E z ) and from a 2.5 
horizontal (E#) 
broadside ideal dipole 
over a perfectly ^ 2 

conducting ground. E z £ 

is a constructively a? 1 

interfered-with field, & 

while Ea, is a S 

destructively o i 

interfered-with field. 
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an antenna that radiated fields polarized such that the electric field was 
essentially perpendicular to the ground plane would produce the largest 
fields and, consequently, the greater power densities on a target or near 
the ground. This is generally not the case for real earth grounds. While 
the perfectly conducting ground plane model may provide some insight 
into the interaction with horizontally polarized fields, it clearly is an inad- 
equate model for vertical polarization. The complex reflection coefficients 
for a plane wave polarized with the electric field in the plane of incidence 
(vertical polarization) and with a plane wave polarized with the electric 
field perpendicular to the plane of incidence (horizontal polarization) are 
given by the following Fresnel expressions. For vertical polarization, it is 
(Collin, 1985) 
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and for horizontal polarization, it is 
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where is the angle of incidence as measured between the ray path and 
the surface of the earth (see fig. 1). These reflection coefficients for the five 
cases of table 1 are presented in figure 3. One can clearly see that for shal- 
low grazing angles, a Brewster angle effect exists (often called a pseudo- 
Brewster angle). The result is destructive interference of the vertically po- 
larized field for incident angles that are less than this pseudo-Brewster 
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Figure 3. Reflection 
coefficients for 
vertically and 
horizontally polarized 
plane waves of 1.3-GHz 
frequency incident on a 
flat ground for five 
earth-parameter cases 
listed in table 1. Note 
that case 5 provides a 
true Brewster angle 
(since conductivity is 
zero) and case 4 is very 
near to producing a true 
Brewster angle at 
1.3-GHz frequency: 

(a) vertical polarization 
and (b) horizontal 
polarization. 





200 
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170 1 1 1 1 

0 20 40 60 80 
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angle. This would lead one to conclude that horizontally and vertically po- 
larized field levels above a ground would be comparable at shallow angles 
of incidence. The Fresnel equations (22) and (23) do not take into account 
the surface or lateral wave produced by finite sources above the ground. 
They are also not applicable to near-field problems. We can compare the 
reflection coefficients as determined from the Fresnel equations with the 
results of our complete solution. To do this with the complete equations 
(eqs (2) to (4) for the vertical dipole and eqs (10) to (15) for the horizontal 
dipole), we fix the radial (or range) distance to 1000 m. The dipole height is 
varied to provide a range of incident angles. The observer will be fixed at a 
1-m height. The incident and "reflected" fields are 


-f'vtotal — \J Ep + E% , 

(25) 

for the vertical dipole and 


-^htotal = \J^4> + Ep + E% , 

(26) 


for the horizontal dipole. For broadside calculations, this simply becomes 
E 0. The terms involving n in equations (2) to (4) and (10) to (15) repre- 
sent the incident field and the remaining terms represent the reflected and 
surface wave terms. The reflection coefficient will be calculated as the to- 
tal incident electric field divided into the total reflected and bound (sur- 
face wave) fields. Figures 4 and 5 provide the comparison of the magni- 
tude of the reflection coefficient for the five cases of table 1. For the most 
part, the Fresnel expressions are a good approximation for the calculation 
of horizontally polarized field values (fig. 5) above realistic earth ground at 
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Figure 4. Fresnel 
reflection coefficient for 
a vertically polarized 
plane wave field 
compared to "reflection 
coefficient" as 
calculated by complete 
formulation for a 
vertical dipole over 
homogeneous earth as 
observed at 1-m height 
1000 m down range at a 
frequency of 1.3 GFIz: 
(a) sea water, (b) wet 
earth, (c) dry earth, 

(d) lake water, and 

(e) dry sand. 
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1.3 GHz. The horizontally polarized fields diverge from the Fresnel expres- 
sions for low conductivities and relative dielectric constants. The resulting 
total fields in such cases predicted by the Fresnel expressions will be larger 
than the complete solution prediction. The results for vertical polarization 
(fig. 4) show significant divergence for shallow incident angles. In this case 
the resulting total fields predicted by the complete solution can be signifi- 
cantly higher than that predicted from employing the Fresnel expressions. 
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Figure 5. Fresnel 
reflection coefficient for 
a horizontally polarized 
plane wave field 
compared to "reflection 
coefficient" as 
calculated by complete 
formulation for a 
horizontal dipole over 
homogeneous earth as 
observed at 1-m height, 
broadside 1000 m down 
range at a frequency of 
1.3 GFIz: (a) sea water, 
(b) wet earth, (c) dry 
earth, (d) lake water, 
and (e) dry sand. 
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Figures 6 and 7 are plots of the fields produced by the two-dipole orien- 
tations; these plots compare the predominant field components for each 
dipole orientation at 1000 and 100 m, respectively. The plots for the 100-m 
case (fig. 7) show the beginning of the lobe effect produced by the phase 
difference between the incident and reflected wave. This is shown more 
dramatically in figure 8. 


Figure 6. Comparison of 
principal fields from an 
ideal dipole oriented 
perpendicular and 
horizontal to a 
homogeneous flat earth. 
In each case, dipole is 
placed 2 m above 
ground plane and 
observer or target is 
1000 m down range: 

(a) sea water, (b) wet 
earth, (c) dry earth, 

(d) lake water, and 

(e) dry sand. 
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Figure 7. Comparison of la > 3 
principal fields from an 
ideal dipole oriented 
perpendicular and S. 2 

horizontal to a H 

homogeneous flat earth. 5 1 5 
In each case, dipole is | 1 

placed 2 m above 
ground plane and 0 5 

observer or target is 0 

100 m down range: 

(a) sea water, (b) wet 
earth, (c) dry earth, 

(d) lake water, and 

(e) dry sand. 2.5 


o> 



0.5 


0 

0 


SIGMA = 4 S/m, EPSREL = 80 
dipole height = 2 m, frequency = 1 .3 GHz 
range = 1 00 m 


Horizontal E-field, horizontal dipole 
z-directed E-field, vertical dipole 


0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Normalized field magnitude (V/m) 


SIGMA = 0.04 S/m, EPSREL = 8 
dipole height = 2 m, frequency = 1 .3 GHz 
= 100 m 


Horizontal E-field, horizontal dipole 
z-directed E-field, vertical dipole 


0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Normalized field magnitude (V/m) 


(b) 


SIGMA = 0.4 S/m, EPSREL = 12 
dipole height = 2 m, frequency = 1 .3 GHz 
2.5 range = 1 00 m 


Horizontal E-field, horizontal dipole 
z-directed E-field, vertical dipole 


J 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Normalized field magnitude (V/m) 


(d) 3 

SIGMA = 0.004 S/m, EPSREL = 80 
dipole height = 2 m, frequency = 1 .3 GHz 
2.5 range = 100 m 


Horizontal E-field, horizontal dipole| 
z-directed E-field, vertical dipole 


0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Normalized field magnitude (V/m) 


(e) 


3 SIGMA = OS/m, EPSREL = 0.2 

dipole height = 2 m, frequency = 1 .3 GHz 
2.5 range = 1 00 m 


2 


1.5 


1 


Horizontal E-field, horizontal dipole 
z-directed E-field, vertical dipole 

0 

0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

Normalized field magnitude (V/m) 


12 





Figure 8. Effect of (a) 

ground reflection on 
primary field 
components near 
ground for typical earth 
parameters. Phase 
difference between 
incident and reflected 
waves results in '© 

SZ 

development of field s 

lobes that are more fe 

(/) 

pronounced as O 

radiating antenna is 
approached: (a) 100 m 
down range and (b) 

1000 m down range. 
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Normalized field magnitude (V/m) 


0.02 
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3. Conclusion 


A suite of MATLAB® m-files have been developed to calculate the electro- 
magnetic fields produced by a vertical and a horizontal inf in itesimal unit 
dipole over a homogeneous flat (ground) plane. Calculations for the fields 
above the ground plane have been made for various ground-plane conduc- 
tivities and relative dielectric constants. The calculations bound the practi- 
cal range of parameters representative of natural earth terrain. 

For a frequency of 1.3 GHz, where the dipole and the observer are close to 
the ground plane (<3 m), significant difference is seen in the magnitude of 
the fields from either dipole orientation. The power density on target will 
be much larger for vertical dipole orientation. The effects of rough terrain, 
foliage, or scattering from manmade or natural objects in the path from the 
dipole to target may alter this conclusion. How these other scatterers might 
affect the field structure at a target at 1.3 GHz is not known at this time. If 
the effects are random in nature, the present conclusion is most likely still 
valid. From the simple, smooth, flat ground model, then, one must con- 
clude that vertical polarization (antenna radiating a vertically polarized 
field) will deliver the most energy to the target. Unless the target's pref- 
erence for field orientation for maximum pickup is known, a vertically po- 
larized antenna may in fact be the best choice for a ground weapon system. 
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Appendix. MATLAB® m-Files 


This appendix documents the MATLAB® m-files that implemented the field 
equations for the vertical and horizontal dipole and array cases. Slight mod- 
ification may be required to obtain all the calculations presented. 


WERF Function 

function w = werf(zj) 

7 WERF(Z,K) Plasma Dispersion Function with complex argument. 

7 

7 Computes the function w(z)=exp(-z"2)*erfc(-iz) using a rational 
7 series with I terms. I should be a power of 2 or it gets SLOW. 

7 Default value of U is 64. z can be a matrix of values. 

7 Taken from Siam Journal on numerical Analysis, Oct 1994, V31 , #5 . 
'/, Modified by R. Chase to work for all z (4/3/95). 

X 

l N=32 gives approx 14 place accuracy, N=64 is better and 
7 it seems as fast. 

X 

7, See wfn.m — Draws graph in Abramowitz k Stegun, 

7 Handbook of Mathematical Functions, p. 298 


if nargin == 1, 1=64 ; end 

7 Default value for U 

11=2*11; 12=2*11; k= [-M+1 : 1 :M-1] ’ ; 

7 M2=no. of sampling points 

L=sqrt (N/sqrt(2) ) ; 

7 Optimal choice of L 

theta=k*pi/M; t=L*tan(theta/2); 

X Variables thetaand t 

f=exp(-t.“2).*(L“2 + t."2); f=[0;f] ; 

7 function to be transformed 

a=real (f ft (f ft shift (f ) ) ) /M2 ; 

7 Coefficients of transform 

a=flipud(a(2:U+l)); 

7 Reorder coefficients 

nz=imag(z) <= 0; 

7 Find im(z) <=0 

z(nz)=conj(z(nz)); 

X Use conj for above 

Z=(L+i*z)./(L-i*z); p=polyval(a,Z); 

7 Polynomial evaluation 

w=2*p./(L-i*z).“2+(l/sqrt(pi))./(L-i*z); 

7 Evaluate w(z) 

w(nz)=2*exp(-(conj(z(nz)).“2)) - conj (w(nz) ) ; 

7 Handle im(z) <= 0 

7,if all(imag(z)==0), w=real(w); end 

7 Rtn real if real 
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Constants 


function [w, cv,epso , uO ,zo ,k0, kl ,kappap ,lo]=cnstdg(f ) 

7 This i-file contains the basic constants to be used by various 
7 m-files and functions 

7 

7 f is input in GHz 
7 global EPSEEL SIGH 
w=2*pi*f*le9; 
cv=2 . 99792458e8; 
epso=(l/(4*pi))*le?*(l/cv"2); 
epsl=EPSREL ; 
u0=4*pi*le-7 ; 
zo=sqrt(uO/epso) ; 
kO=w*sqrt (epso*uO) ; 
epsr=(epsl-j*SIGMA./(»*epso)); 
kl=kO*sqrt(epsr) ; 
kappap=-kO/sqrt (epsr+1 ) ; 
lo=cv/(f *le9) ; 


Vertical Dipole 

function [ez ,er,hp]=dipg(f ,d,rho,z) 

7 

7 Electromagnetic fields from a vertical z-directed current element 
7 at a height, d, over a conducting dielectric plane (ground) calculated 
7 by the King/Sandler model This formulation is that developed 
X as equations 6, 7, and 8 of "The Electromagnetic Field of a Vertical Electric 
7 Dipole over the Earth or Sea", IEEE Trans, on Antennas and Propagation, 

7 March 1994, Vol 42 Bo. 3, page 383. It is the same as the King/Owens/Wu 
7 formulation. This formulation is that developed as equations 4.2.30-32 of 
7 "Lateral Electromagnetic Haves", Springer-Verlag , 1992, pp 293-297. 

7 

X > f- frequency; rho- radial distance from the z-axis 

7, < z- height of observer above ground plane (input as an array) 

X > ez,er,hp- field components at the observer, where the second letter 

7 < designates the component — z, r (rho), p (phi) 

7 

7 The formulation is subject to Ikl(ground) |>3|k2(air) | 

7 

[w, c ,eps0 ,u0 ,zo ,k2,kl ,kappap, lo] =cnstdg(f ) ; 

X 
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7 King, et. al assume an exp(-i*omega*time)dependency thus we convert kl 

7 

kl-conj(kl); 

7,JI2=(kl/k2)"2; 
k21=k2/kl ; 

7 

7 He assume a unit dipole 

7 

il=l ; 

rl=(rho“2+(z-d) . '2) . " . 5 ; 
r2=(rho“2 + (z + d) . "2) . “ . 5 ; 
sd=rho . /rl ; 
sr=rho7r2; 
cd=(z-d)7rl; 
cr=(d+z)7r2; 

p z (r2*k2“3/(2*kl“2))^((k2*r2+kl*(z+d))7(i2*rho))/2; 

7 

7 The attenuation function less the exp(i*p) term - related to the 
7 error function is 

7 

Ferf=((l+i)/4)*werf (sqrt(i*p)) ; 

7 

7 Develop the terms for the field expressions 

7 

tl=w*u0*il/(2*pi*k2); 
tll=-il/(2*pi); 
t2=exp(i*k2*rl )/2 ; 

t3 = (i*k27rl)-(t7rl."2)-i7(k2*rl.~3); 

t31=(i*k27sfi-(l7rl.“2); 

t4=(cd.''2).*((i*k27rl)-(37rl.‘2)-(3*i) 7(k2*rl.'3)) ; 

Ml=((i*k27rl)-(3. /rl . ~2)-(3*i) . /(k2*rl . "3)) ; 
t5=exp(i*k2*r2)/2; 

tMi*k27r2)-(l7r2."2H7(k2*r273); 

t61=(i*k27r2)-(L/r2.“2); 

t7=(cr.‘2).*((i*k27r2)-(37r2.~2)-(3*i) 7(k2*r2.~3)); 
t71=((i*k27r2)-(37r2.'2)-(3*i) 7(k2*r2. ‘3)); 
t8=(2*t5*(k2"3)/M') . *((pi.-/(k2*r2)). “ . 5) . *sr. *Ferf ; 
t8M(k2"3)/kl)*((pi7(k2*r2)).".5).*Ferf; 
t82=(2*t5 . *(k2~3)/kl) . *( (pi 7(k2*r2)) . " . 5) . *Ferf ; 

7 

7 Calculate the fields 
7 ez=tl*(t2. *(t3-t4)+t5. *(t6-t7)-t8) ; 
er=-tl*(t2.*sd>cd.*t41+t5.*sr.*cr.*t71-k21*2*t5.*(sr.*t61-t81)); 
hp=tll*(t2.*sd.*t31+t5.*sr.*t61-t82); 
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Horizontal Dipole 


function [ez ,ep,er J hz,hp,hr]=dipgh(f ,d,phi ,rho ,z) 

'/, Electromagnetic fields from a horizontally directed (perpendicular to z) 

'/, current element located in the phi=0 plane 

'/, at a height, d, over a conducting dielectric plane (ground) calculated 
'/, by the King/Owens/Hu formulation. This formulation is that developed 
•/, as equations 7.10.75, 7.10.80, 7.10.84, 7.10.92, 7.10.93 and 7.10.94 of 
'/, "Lateral Electromagnetic Waves" , Springer-Verlag , 1992, pp 293-297. 

'/, > f- frequency(GHz); phi- angle about the z-axis, rho- radial distance from 

'/, < the z-axis, z- height of observer above ground plane (input as an array) 

'/, > ez,ep,er,hz,hp,hr- field components at the observer, where the second 

'/, < letter designates the component — z, r (rho), p (phi) 

X 

'/, 

'/, The formulation is subject to Ikl(ground) |>3|k2(air) | 

[w, c ,eps0 ,u0 ,zo,k2, kl .kappap, lo] =cnstdg(f ) ; 

'/, King, et. al assume an exp(-i*omega*time)dependency thus we convert kl 

kl=conj (kl) ; 

, /,I2=(kl/k2)"2; 
k21=k2/kl ; 
k2p=k2/rho; 

'/, Se assume a unit dipole 

’/, 

il=l ; 

cp=cos(phi) ; 
sp=sin(phi) ; 

rl=(rho“2+(z-d).“2).“.5; 

r2=(rho~2+(z+d).~2).~.5; 

sd=rl/rho ; 

sdo=l . /sd; 

sr=r2/rho ; 

sro-l./sr; 

cd=(z-d) . /rl ; 

cr=(d+z)./r2; 

p-(r2*k2“3/(2*kr2)).*((k2*r2+kl*(z+d))./(k2*rho)).“2; 
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'/, The attenuation function less the exp(i*p) ten - related to the 

'/, error function is 

'/, Ferf=((l+i)/4)*werf (sqrt(i*p)) ; 

'/, 

'/, Develop the terms for the field expressions 

'/, 

tl=w*u0*il*cp/(4*pi*k2) ; 
tll=-w*u0*il*sp/(4*pi*k2) ; 
tl2=il*sp/(4*pi) ; 
tl3=il*t|/(4*pi); 
t2=exp(i*k2*rl); 

t32=(i*k2./rl)-(l./rl."2)-i./(k2*rl."3); 
t3°(2. /rl . ~2)+2*i 7(k2*rl . “3) ; 

t31=(i*k2./rl)-(l./rl.“2); 

H=(cd.^).*((i*k2./rl)-(3./rl.“2)-^i)./(k2*rl.“3)); 

t41=(cd).*((i*k2./rl)-(3./rl.“2)-(3*i)./(k2*rl."3)); 

t5=exp(i*k2*r2); 

t62=(i*k2./r53-(l./r2.~2)-i./(k2*r2."3); 

t6=(2./r2.“2)+2*i./(k2*r2.“3j; 

t61=(i*k2./r2)-(l./r2.“2); 

t7=(cr.‘2).*((i*k2./r2)-(3./r2.‘2)-(3*i) ,/(k2*r2.“3)) ; 
t71=(cr).*((i*k2./r2)-(3./r2.“2)-(3*i)./(k2*r2.'‘3)); 
t72=((l-/r2.~2) + (3*i)./(k2*r2."3)-3./((k2~2)*r2.~4)); 
t73=(cr.“2).*((i*k2./r2)-(6./r2.“2)-(15*i)./(k2*r2.“3)); 
t8=((k2“3)/kl)*((pi . /(k2*r2)).“.5) . *sr.*Ferf ; 
t81=((k2“3)/kl )*( (pi . /(k2*r2) ) . “ . 5) . *Ferf ; 
t82=(2*t5*(k2“3)/kl)>.(fpi./(k2*r2)).“.5).*Ferf; 
t83= ((pi . /(k2*r2) ) . ~ . 5) . *Ferf ; 

•/. 

I Calculate the fields 

’/, 

ez=tl*(-t2.*sdo.*t41+t5.*sro.*t71-2*k21*t5.*(sro.*t61-t81)); 
er=t l*(t2 . * (t3+t4)-t5. *(t6 + t7-2*k21*t61+2*(k21 ~2)*(t62-t,8) ) ) ; 
ep=tll*(t2.*t32-t5.*t62-t5.*(-2*k21*cr.*t61+2*(k2r2)*(t6 + t7). . . 
+2*i*k2r3*k2p*(sr.~2).*t83)); 

hr=tl2*(t,2. *cd.*t31-t5.*cr. *t61+2*k21*t5 . *(t6+i*k21*k2p*(sr . ‘2) . *t83+t7)) 
hp = tl3*(t2.*cd.*t31-t5.*cr.*t61 + 2*k21*t5.*(t62-k21*(k2“2)*(sr).*t83)); 
hz=tl2*(t2. *sdo. *t31-t5. *sro. *t,61+2*sro. *t5.*(k21*t71(k21'2)*(t72+t73))) ; 



Plot m-File for Fields 

'/, 

'/, This m-file plots the fields over a conductive flat earth produced by an ideal 
'/, dipole placed a distance d above the earth. It compares the results from 
'/, a vertical and horizontal dipole. 

'/, Establish the problem conditions 

'/, EPSREL- Relative dielectric constant; SIGMA- Earth conductivity (S/m) 

EPSREL=80 ; SIGMA=4 ; 
global EPSREL SIGMA 
, /,EPSREL= 1 2 ; SIGMA= . 4 ; 

'/,EPSREL=8 ; SIGMA= . 04 ; 

'/,EPSREL=80 ; SIGMA= . 004 ; 

‘/,EPSREL=2 ; SIGMA= . 000 ; 

'/, Location of dipole (m) 

d=2; 

'/, Location of observer, rho (m); phi (radians); z (m) — > an array 

rho-1000; 
phi=pi/2 ; 
z=. 1*[1: 1:30] ; 

'/, f- frequency in GHz 

f=l .3; 

'/, Field normalization factor - one Watt radiated 

'/, 

[w, c. ,eps0 ,u0 ,zo ,k2, kl ,kappap, lo] =cnstdg(f ) ; 

Fn=sqrt(12*pi/((k2“2)*zo)); 

'/, Horizontal dipole fields 

[ez , ep , er , hz , hp , hr] =dipgh (f , d , phi , rho , z) ; 

plot(Fn*abs(ep),z,'-b’) 

hold 
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'/, Vertical dipole fields 

'/, 

[ezv,erv,hpv]=dipg(f ,d,rho,z); 
plot (Fn*abs(ezv) , z, ' - . r ; ) 
ts=Fn*abs(erv(l)); 

title( 'Vertical and Horizontal Ideal Dipole fields over Ground’) 
text(ts,z(19), [’SIGMA - ’ ,mnn2str (SIGMA) , ’ S/m EPSREL- 
’ ,num2str(EPSREL)]) 

text (ts,z(18) , [’Dipole height = ’ ,nran2str(d), ’ meters; Frequency z 
’ ,num2str(f), ’ GHz’]) 

text (ts,z(17) , [’Range - ’ ,num2str (rho) , ’ meters’]) 

ylabel ( ’ Observer height - meters’) 

xlabel( ’formalized field magnitude - V/m’ ) 

legend(’Horizontal E-field, horizontal dipole ’ , ’z-directed E-field, 

vertical dipole’) 

hold 

Fresnel Reflection Coefficients and Plots 

function [rcomv,rcomh]-Fresnell(f ,psi) 

•/, 

'/, This routine calculates the Fresnel reflection coefficients for 
'/, vertical and horizontal incident fields 

'/, 

'/, after Collin, "Antenna and Radiowave Propagation", page 345 

'/, f-frequency in GHz; psi- array of incident angles 

'/, Code returns the arrays rcomv and rcomh, the vertical and horizontal 
'/, complex reflection coefficients, respectively. 

’/, 

global EPSREL SIGMA 
w=2*pi*f*le9; 
cv=2 . 99792458e8 ; 
epso=(l/(4*pi))*le7*(l/cv~2); 
epsr=(EPSREL-j +SIGMA . / (w*epso) ) ; 

rcomv z (epsr*sm(psi)-sqrt (epsr-cos(psi) . “2)) . / (epsr*sin(psi)+sqrt (epsr-cos(psi) . “2)) ; 
rcomh = (sin(psi)-sqrt(epsr-cos(psi).“2))./(sin(psi)+sqrt(epsr-cos(psi).“2)); 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



'/, This m-file plots the reflection coefficient for vertical and horizontal 
'/, fields over a homogeneous ground plane 

'/, EPSREL- Relative dielectric constant; SIGMA- Earth conductivity (S/m) 

EPSREL=80 ; SIGMA=4 ; 
f*1.3; 

global EPSREL SIGMA 
psi= [pi/1000 :pi/1000:pi/2] ; 

[rcomv,rcomh]=Fresnell(f ,psi) ; 

EPSREL=12 ; SIGMA= . 4 ; 

[rcomvl ,rcomhl]=Fresnell (f ,psi ) ; 

EPSREL=8 ; SIGMA= . 04 ; 

[rcomv2,rcomh2]=Fresnell (f ,psi ) ; 

EPSREL=80 ; SIGMA= . 004 ; 

[rcomv3,rcomh3]=Fresnell (f ,psi); 

EPSREL=2 ; SIGMA= . 000 ; 

[rcomv4,rcomh4]-Fresnell(f,psi); 

subplot (2 ,2 , 1) .plot ( (180/pi)*psi , abs(rcomv) , J k 1 ) 
hold 

axis( [0, 90 ,0,1]) 

xlabel(’Incidence angle (degrees)’) 
ylabel( ’Reflection Coefficient’) 

'/,text(5, . 8 , [’SIGMA = ’ ,num2str(SIGMA) , ’ S/m EPSREL- ’,num2str(EPSREL)]) 

subplot (2,2, 1) .plot ( (180/pi)*psi , abs(rcomvl), ’r’) 

subplot (2,2,1), plot ( (180/pi)*psi, abs(rcomv2), ’b’) 

subplot (2, 2,1), plot ( (180/pi)*psi , abs(rcomv3), ’g’) 

subplot (2, 2,1), plot ((180/pi)*psi,abs(rcomv4), ’c’) 

hold 

subplot (2 ,2,2) ,plot ( (180/pi)*psi , (180/pi )*angle(rcomv) , ’k’ ) 
hold 

axis( [0,90,-190,5]) 

xlabel ( ’ Incidence angle (degrees)’) 

ylabel( ’ Phase ’ ) 

subplot (2 ,2 , 2) .plot ( (180/pi)H<psi , (180/pi )*angle(rcomvl) , ’r ’ ) 
subplot (2 ,2 , 2) ,plot ( (180/pi)*psi , (180/pi)*angle(rcomv2) , ’b ’ ) 
subplot (2 ,2 , 2) ,plot ( (180/pi )*psi , (180/pi)*angle(rc:omv3) , ’g ’ ) 
subplot (2, 2,2), plot ( (180/pi)*psi , -(180/pi )*angle (rcomv4) , ’ c ’ ) 
hold 

subplot (2,2,3), plot ( (180/pi)*psi , abs(rcomh) , ’k ’ ) 


24 



hold 

axis( [0,90, .5,1]) 

xlabel ( 1 Incidence angle (degrees)’) 

ylabel( ’Reflection Coefficient’) 

subplot (2,2,3), plot ( (180/pi)*psi, abs(rcomhl), ’r’) 

subplot (2,2,3), plot ( (180/pi)*psi, abs(rcomh2) , ’b’) 

subplot(2,2,3) ,plot ( (180/pi )*psi ,abs(rcomh3) , ’g') 

subplot (2 ,2 ,3) ,plot ( (180/pi)*psi , abs(rcomM) , ’ c ’ ) 

hold 

subplot (2, 2, 4) .plot ((180/pi)*psi , (180/pi)*angle(rcomh) , ’k’ ) 
hold 

axis( [0,90,170,200] ) 
xlabel(’Incidence angle (degrees)’) 
ylabel ( ’Phase ’ ) 

subplot (2 ,2 ,4) .plot ( (180/pi )*psi , (180/pi) *angle(rcomhl ) , ’r ’ ) 
subplot (2 ,2,4) .plot ( (180/pi)*psi , (180/pi )*angle(rcomh2) , ’b ’ ) 
subplot(2,2,4),plot((180/pi)*psi,(180/pi)*angle(rcomh3), ’g’) 
subplot (2 , 2 ,4) .plot ( (180/pi)*psi , (180/pi )*angle(rc.omh4) , ’ c ’ ) 
hold 
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